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A detailed quantum kinetic master equation is developed which couples the kinetics of a trapped conden- 
sate to the vapor of non-condensed particles. This generalizes previous workwhich treated the vapor as being 
undepleted. 
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I. INTRODUCTION 

Until the experimental realization of a Bose-Einstein con- 
densate of magnetically trapped Alkali atoms [[l|-^ in 1995 
(which has now become a widely available technology [^) 
the theory of Bose-Einstein condensation was dominated by 
the desire to understand the behavior of superfluid liquid he- 
lium, of which it was assumed that the weakly interacting 
degenerate Bose gas was a simplified and admittedly inade- 
quate model. The central tools for the description of BEC 
were the Bogoliubov theory of superfluidity and modifications 
of it [H"^, the Gross-Pitaevskii equation [|^] and techniques 
based on the Green's function formalism of many body the- 
ory [ pp. T his philosophy is excellently presented in Griffin's 
book[0. 

The prospect of actually producing a degenerate Bose gas 
gave rise to a number of investigations [|l9|| based largely on 
a combination of kinetic theory as described by the Uehling- 
Uhlenbeck equation [|l8|] and the the Gross-Pitaevskii equa- 
tion. However, nearly all of these investigations were from the 
point of view of macroscopic statistical mechanics, in which 
the BEC was seen as a sample of an extended fluid, rather 
than the tightly trapped condensate, which behaves more like 
a very large atom or molecule than a droplet of macroscopic 
fluid. 

Our own program of quantum kinetic theory has been for- 
mulated in such a way as to take advantage the simplifications 
which can be made for a tightly trapped condensate. It has 
been partially presented in our first four papers on the subject 
|[ri|], as well as in two papers on the growth of a Bose-Einstein 
condensates [ll6(T7|. However, all of the previous papers cov- 



ered only special situations. In QKI we treated only spatially 
homogeneous systems with at most a small amount of con- 
densate; QKII dealt with stochastic equations as a model of 
the initiation of the condensate, QKIII and QKIV considered 
only the case where the majority of the atoms provided a bath 
of atoms with energy above a certain value (called En) of 
fixed chemical potential fi and temperature T for a relatively 
small number of lower energy atoms, called the condensate 
band, which included the condesnsate level itself as well as 
the levels with energy less than En. 

In this paper we will present what we see as a compre- 
hensive and reasonably practical description of the subject, 
in which the various restrictions in the previous papers are 



largely eliminated. The description which results is best de- 
scribed in terms of a system composed of a condensate and a 
vapor in interaction with each other. The vapor is described 
by a quantum kinetic master equation of the form introduced 
in QKI, but appropriately adapted to take account of the pres- 
ence of a trap. For nearly all situations, this is equivalent to a 
quantum Boltzmann equation of the Uehling-Uhlenbeck | ]l8[ ] 
form. The condensate is described by a fully quantum me- 
chanical Hamiltonian, which for many purposes can be treated 
by a Bogoliubov approximation, but truncated to include only 
those excitations not already included in the vapor. The in- 
teraction between vapor and condensate consists of two parts; 
a Hamiltonian mean field term, and a master equation term 
which describes the transfer of energy and particles between 
vapor and condensate. The largest part of the mean field term 
can be estimated, and included by defining effective poten- 
tials for the condensate and the vapor. The computation of 
these effective potentials involves a procedure very similar 
to that of the Hartree-Fock-Bogoliubov [|2^] method. How- 
ever these effective potentials are used only to give a basis 
set of eigenfunctions, in terms of which we can describe the 
kinetics; the resulting kinetic equations give rise to correc- 
tions to the corresponding eigenvalues — our use of a basis 
of eigenfunctions similar to those given by the Hartree-Fock- 
Bogoliubov method does not mean our method is the same as 
the that method, since there is a residual term remaining after 
the majority of the mean field term has been included in the 
effective potentials, and this is of a similar order of magnitude 
to the master equation term. 

This paper is concerned mainly with methodology and rea- 
sonably careful derivations, rather than with directly applica- 
ble results. However, it provides a justification for the heuris- 
tic description of the theory of condensate growth presented 
in [[l7|], and does provide a way in which hydrodynamics and 
quantum kinetics can be connected together. There are, of 
course alternative approaches to the general issue of combin- 
ing kinetic theory and quantum mechanics, which have dif- 
ferent emphases. For example Stoof [|T]| has chosen to de- 
velop an approach which avoids the use of the pseudopoten- 
tial method, which he has shown will fail very near the criti- 
cal point; Walser et al. [ p2| ] have developed a kinetic approach 
somewaht more related to ours; while Zaremba et al. \ 23 1 have 



developed a two fluid hydrodynamic description. The major 
criterion for choosing a particular approach will ultimately be 
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the degree of usefulness in explaining experiments, and its 
ease of use — but we will leave this judgment on this to the 
BEC community. 

In the remainder of this introduction we will summarize the 
philosophy and scope of our methodology, and to emphasize 
the most significant results. 

A. The field theory of Bose particles 

We want to consider the Bose atoms to be described by a 
second-quantized field, in the pseudopotential approximation; 
that is, we write 

H = H^,,,+Hi + Ht, (1.1) 

where 

Hi = ^Jd^xJ d3x'V't(x)?/.t(x')M(x-x')?A(x')^(x). 

(1.3) 

and the term Ht arises from a trapping potential as 

iJr = y"rf^xV"T(x)^A^(x)V'(x). (1.4) 

In this paper we will not go deeply into the approximations 
which allow us to use the pseudopotential form of the Hamil- 
tonian, in which 

Zl 77-/7 

w(x' - x) ^ ii(5(x' - x) = (5(x'-x), (1.5) 

m 

which is thoroughly investigated in many works [p4||. Suffice 
it to say only that this is an approximation whose validity is 
assured provided the wavelengths of the atoms are very much 
longer than the range of the potential u(x' — x), and provided 
that we are not working too close to the critical point. Thus, 
at typical temperatures and densities currently used, this ap- 
proximation is valid either when there is a significant amount 
of condensate, or when there is no condensate. The essence 
of the derivation of this pseudopotential approximation is re- 
ally an adiabatic elimination of the fast motion of the atomic 
wavefunctions during the very small proportion of time during 
which they are actually undergoing collisions. Thus the result 
is effectively a Hamiltonian which is coarse-grained in time, 
and thus, because of the dispersion relation, also in space. The 
details of this philosophy are available in a number of refer- 
ences [ ^4] ] Having made this approximation, we know we are 
then justified in using an approach which is essentially pertur- 
bative. 



B. Thermalization assumptions 

It is evident in experiments that a typical condensate in a 
trap is accompanied by a thermal vapor cloud. The conden- 
sate consists of a large number of particles, all of which share 



the same wavefunction, while the vapor is essentially com- 
pletely thermalized — it thus consists of particles which do not 
share wavefunctions with each other, and indeed most wave- 
functions are not occupied by any particles at all. The va- 
por consists mainly of the higher energy particles, but there 
is a transition from vapor to condensate which, althought not 
sharp, necessitates the drawing of a definited boundary be- 
tween them. 




Fig. 1: Representation of the modification of the trapping potential 
for a) the non-condensate band and b) the condensate band due to 
mean-field effects 



1. Condensate and non-condensate bands in the previous papers 

The concept of condensate and non-condensate bands was 
introduced by us in QKllland [|l|]. In these earlier works we 
chose to make the distinction between these in terms of the ef- 
fect of the presence of a condensate on the trap energy levels. 
This meant that the non-condensate band was taken to con- 
sist of those energy levels above a certain value, E^t, which 
were treated as being essentially unaffected by the conden- 
sate; the condensate band consisted of all levels with energies 
E < Efj, and these were taken to have wavefunctions and 
energies which were given by the Bogoliubov method. 

In this paper we move to a more useful way of treating 
the two bands. In QKI we developed a wavelet descrip- 
tion of weakly condensed systems, which allowed for a fully 
quantum-mechanical description of the thermalization, but 
this was derived only for the case of free particles. It was 
the absence of a good wavelet description of trapped particles 
that led us to consider in QKIII only the approximation that 
the non-condensate band could be treated as being like a bath 
of thermal particles with a fixed chemical potential, /i and a 
fixed temperature, T. Attention was then concentrated on the 
derivation QKIII and successful application [ p^p3| , p^ of a 
master equation for the condensate band. 
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2. Condensate and non-condensate bands in this paper 

In this paper we have worked out how the wavelet descrip- 
tion can be modified in the presence of a trapping potential, 
and have been able to give a full set of coupled equations con- 
necting the two bands. The distinction between the two bands 
is now made in terms of whether the energy level is particle 
like, or phoiion-like. We take the division between the bands 
to be at En, where Eb, is the energy above which all the exci- 
tations can be considered, to a good degree of approximation, 
to consist of a definite number of particles, unlike the phonon- 
like excitations, which have lower energies. 







Nj, increasing 



large 



Fig. 2: Representation of the non-condensate-band (left) and con- 
densate band (right) effective potentials appropriate to an increasing 
number of condensate atoms. Fore each particular sector, defined by 
the numbers A^^nc and Nc, tiie different effective potentials define 
the basis eigenstates in terms of which the master equation is de- 
fined. The process of growth is that of transition from one sector to 
the next, and does not involve a time dependence of the potential. 

In order to take account of mean-field effects of the 
condensate-band atoms on the non-condensate-band atoms, 
and conversely, we use effective trap potentials, for the con- 
densate and non-condensate bands which include these mean- 
field effects. Thus, we write 



l/j^c'W = 2u[p(x)+fi(x) 
V^\x) = 2un(x), 



(1.6) 
(1.7) 



where p(x) is the mean condensate band density, and n(x) 
is the mean density of the vapor in the non-condensate band. 
The values of these quantities depend on the particular phys- 
ical situation, including the numbers of particles, A^nc ™d 
Nc, in the respective bands. 
The effective potentials are then 



yc^«(x) 



^t(x) 
^t(x) 



^NC(X) 



(1.8) 
(1.9) 



As written, the effective potentials contain no time- 
dependence. However, time-dependent situations can be dealt 
with in two ways. 

Firstly, consider condensate growth. The picture put 
forward by us in [ p"6|p^ ] treats condensate growth as the 
slow change of the populations of the condensate and non- 
condensate bands by transitions of particles from one band 
to the other The transitions are between energy levels deter- 
mined by the values of the densities n(x), p(x) appropriate 
to the actual numbers TYnc, of non-condensate-band and 
condensate-band atoms. Thus the values of ?l,(x) and p(x) 
(and hence also the mean-field ptentials) will change adiabat- 
ically as the condensate grows, and the energy levels needed 
to compute the transition rates will also change adiabatically. 
The basic requirement is that the growth is sufficiently slow 
compared to the enrgy thermalization time within each band 
that there are always well-defined values of ?l(x) and p(x). 

The second way in which time-dependence can be treated 
is perturbatively. If the condensate and vapor oscillate with 
small amplitude about mean values, we can show that this 
can be taken into account by using appropriate time depen- 
dent values of n(x) and p(x). 

These two effective potentials provide a basis within which 
we can develop the full quantum kinetic theory. This is done 
by dividing the field operator into two parts 



ipix) = 0(x) + ■i/'Nc(x) 



(1.10) 



which then correspond to the condensate band and non- 
condensate band degrees of freedom. We then write the 
Hamiltonian ( 1 . 1 ) as 



H = H 



off 



(H - H-*") 



(1.11) 



in which H'^^ is the sum of the condensate-band and non- 
condensate band effective Hamiltonians; hence the non- 
condensate-band effective Hamiltonian is 



frcft 



2m 



y^g(x)U(x) 



(1.12) 



and the condensate band effective Hamiltonian is 



d3x0t(x) 



2m 



+ yc^^(x) 



+ - / d3x0t(x)0t(x)0(x)0(x). 



(1.13) 



The quantities p(x) and n{x.) are nominally the mean 
condensate-band and non-condensate-band densities, but the 
real criterion for their choice is to minimize the effective in- 
teraction term H — 
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3. Basis states for the master equation 

The non-condensate-band effective Hamiltonian is now 
written in terms of a wavelet basis, as follows. We write 



V'Nc(x) = ^ WQ(x,r)AQ,r ^ V^q(x). 
Q,r Q 



(1.14) 



Here Aq ,. are destruction operators with the commutation re- 
lations 



(1.15) 



so that r and Q are discrete indices. The wavelets wq(x, r) 
are a complete orthonormal set, such that 

1. wq(x, r) ^ when |r — x| is large. 

2. wq{x, r) is a linear combination of trap eigenfunctions for 
the effective potential Vt{x) + V^Q{r) with energy eigenval- 
ues in a range restricted to a neig hborhood of Ti^Q^/Zm. 

3. With each set of Q, r we can associate a momentum 
?iK(Q, r) which is such that 



2m 



2m 



(1.16) 



For every direction of Q we can also associate a direction of 
K; for example in a spherically symmetric harmonic trap Q 
and K are parallel. 

In general a rule of association between the directions is not 
simple, but it can be specified implicitly. 

4. Using the function K(Q, r) we can write the commutator 

[V'q(x)>Jj,(x')] « <5QQ.e^K(Q'^)-(---')5Q(x,x') 

(1.17) 

where the function gQ(x, x') is only significantly different 
from zero when x w x', and is normalized to one in both 
variables x and x'. 

This set of wavelets can be explicitly constructed in a num- 
ber of cases, and we assume its existence in general. It pro- 
vides a resolution of the non-condensate field into discrete 
"cells" in phases space labeled by Q and r. The wavelet 
function WQ.r(x) corresponds to a particle approximately lo- 
calized at r in space and with momentum in a small range 
around ?iK(Q, r). Its total energy is thus approximately 
E{Q) = ?i^Q^/2m, which is independent of r. A particle 
in a wavelet state Q, r will evolve under the influence of the 
Hamiltonian -ff^c ^ linear combination of wavelets with 
the same Q, so that the evolution represents energy conserv- 
ing and particle conserving flow in phase space. We thus ob- 
tain a Liouvillian flow in this discrete phase space. 

The evolution of the condensate band is described by 



(1.13), which differs from the full Hamiltonian only in that 
1. The mean-field potential of the vapor Vq^^{x.) occurs. 



2. The operator (f>{x) is limited to the lower energy levels of 
the full Hamiltonian, and thus 0(x) has a non-local commuta- 
tion relation, the degree of non-locality being of the same size 
as the average spatial width of the wavelets used to construct 
the non-condensate-band. 

Although we have described the wavelets by approximate 
properties, once the wavelets have been constructed, the de- 
scription in terms of them is exact. However, in order to use 
this description some approximations are still necessary. The 
simplest procedure — and the one we shall use here — is to use 
the number-conserving Bogoliubov model, but include the va- 
por mean-field potential V™^{x.) ~ 2ufi(x) as well as the trap 
potential. This provides us with an approximate resolution of 
the non-condensate-band field operator in the form 

(/)(x) = A (^^^Jx) + ^_l_ x(x)j (1.18) 

in which 

1 . A is a destruction operator for the total number of atoms in 
the condensate band; 

2. Nc = Av& the number of atoms in the condensate band; 

3. i^Nc (x) is the condensate wavefunction. 

4. The quantity x(x) is a quasiparticle field operator, which 
has the expansion 

X(x) = {p,„(x)5„ + g„(x)6j„} . (1.19) 

TCI 

The (x) and (x) are quasiparticle amplitudes for the 
quasiparticle creation and destruction operators The 
condensate has energy ^(iVc), and each quasiparticles carries 
an energy hem{Nc). 

We assume that Nq is almost equal to the number of parti- 
cles N in the condensate level. This is a significantly weaker 
requirement than that which would be necessary to use the 
usual Bogoliubov method for the whole system, which would 
require the total number of particles in the system, Nq + A^nc > 
to be close to the number of particles in the condensate. 

At this degree of approximation, the description is rather 
like the conventional Hartree-Fock-Bogoliubov |20| method. 
However, there is so far no thermalization assumed — this 
can only be done satisfactorily by introducing irreversibility, 
which we shall do using our master equation. 



4. The master equation 



The master equation has two main purposes 

1 . To take account of the interactions that are included in the 
term H — H"^^ in the equation ( 1.12 ) 



2. To introduce irreversibility, and thus thermalization. 

The technical methodology we use is very similar to that 
used in QKl. The equation of motion for the density opera- 
tor is developed perturbatively by projecting onto the space 
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of density operators in which the relative phases of non- 
condensate-band states with different Q are eHminated. We 
cannot eHminate the phases corresponding to the same Q and 
different r, since the Liouvillian transport generated by is 
off-diagonal in r, would itself be eliminated by such a proce- 
dure. But Liouvillian transport generates the streaming terms 
in the quantum Boltzmann equation, and must therefore be 
preserved if we are to obtain a full hydrodynamic description 
of the vapor 

As in QKI and QKIII the derivation of the master equa- 
tion relies on a perturbation expansion in the interaction term, 
and on a random phase approximation in certain of the coef- 
ficients. The master equation which results is natural, self- 
consistent and consistent with statistical mechanics. 

For example, a stationary solution for the density operator 
can be written 



exp 



h: 



cff 



kT 



exp 



kT 



Here 



h: 



cff 



(1.20) 



(1.21) 



and is given by ( 1.13 ). 

There are thus no correlations between the condensate and 
non-condensate bands, and none between the different Q, r 
"cells" in the non-condensate-band. This description is self- 
consistent if the highest energy states in the condensate-band 
(and all higher energy states) are essentially particle-like, 
since the density operator will then look much the same at 
the top of the condensate-band as at the bottom of the non- 
condensate band. The master equation in the full form is writ- 
ten in Sect. ^IIE , and comprises terms which account for all 
physical processes as follows. 



1. Condensate-band term: This is given by ( |3.35[|3.37[), and 
arises from the condensate band Hamiltonian H^, (1.13). 

2. Non-condensate-band terms: These are given by (3.38a- 



3.38d), and are very like those in QKI. They account for Li- 



ouvillian flow and collisions within the non-condensate-band. 

3. Mean-field coupling term: One part of the interaction 
leads to a Hamiltonian coupling between the two bands, and 
this is a residual part of the Hamiltonian, which can be written 



Q 



(1.22) 



Since this is only a residual after subtracting the mean-field 
terms, it is expected to be of the same order of magnitude as 
the irreversible terms. It is neglected entire ly in writing the 
stationary density operator in the form ( 1.20 ). 

4. Irreversible coupling terms: These are the irreversible 
terms which give rise to energy and particle transfer between 
the condensate-band and the non-condensate-band. They are 
written in terms of some formal condensate-band operators, 
and the field operators ipQ (x) of the non-condensate-band. In 
order to use the master equation, the eigenf unctions of H^, 
the effective condensate-band Hamiltonian, must be found. 
This is a question of choosing the appropriate approximation 
method. 



C. Practical methods 

By assuming that the non-condensate band can be main- 
tained in a thermal stationary state we recover the method- 
ology of QKIII. The next obvious step is to take a time- 
dependent local equilibrium approximation for the non- 
condensate-band, and a number-conserving Bogoliubov ap- 
proximation for the condensate-band. We can then talk about 
the system as a thermalized vapor in the non-condensate-band, 
in interaction with a system composed of a condensate and 
quasiparticles in the condensate-band. This is the most useful 
way of expressing our results. We find that 

1 . The condensate-band is described as the Bogoliubov sys- 
tem of quasiparticles and condensate, in which the mean-field 
potential induced by the vapor may be time-dependent. 

2. The vapor is described by a quantum Boltzmann equation, 
in which the mean field potential induced by the condensate 
band may be time-dependent. 

3. There are transfer terms between the condensate-band and 
vapor, which may be evaluated in various degrees of approxi- 
mation. 



1. Non-condensate-band equations of motion 

The equation of motion for the phase-space density ./k(x) 
of the non-condensate-band vapor can be written 



a/K(x) _ fhK-w^ w^v^U^,t) ■ Vk 



dt 



\ m 

|2 



/k(x) 



+^ / d""^^ J ^'^3 J d3K45(K + K2 - K3 - K4)<5(w + a;2 - W3 - W4) 

Rnc Rnc Rnc 

x{/k(x)/k.(x)[/k3(x) + 1][/k,(x) + 1] - [/k(x) + 1][/k.(x) + 1]/k3(x)/k,(x)} 
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9/k(x) 

at 



+ 



g/K(x) 

dt 



9/k(x) 
dt 



(1.23) 



The final three terms represent transfers of energy and atoms 
to and from the condensate band in which respectively 1, 2 
and 3 condensate-band field operators are involved in the ma- 
trix elements. 



The transfer terms are given in full in (6.12), (6.22) and 
( |6.25 ), and are necessarily rather complicated because of the 
mixing of creation and destruction operators generated by the 
Bogoliubov method. However, the essence of the result is 
that there are simple rate processes involving non-condensate- 
band particle numbers, condensate band quasiparticle num- 
bers, and the number of particles in the condensate. All of the 
rates involve integrals like 

Jd^K2 jd^Ka jd^kS {L02Ki{^) - uj) S{K + K4 - K2 - k) 
X [1 + ./k(x)] [1 + /k,(x)] /K.(x)W/(x,k). (1.24) 



The quantity W/ (x, k) is a Wigner function of a quasipar- 
ticle amplitude — we see therefore that there are energy and 
momentum conservation delta functions as in the quantum 
Boltzmann equation, and where a particle is in the condensate 
band we replace its phase space density by the Wigner func- 
tion corresponding to the particular quasiparticle amplitude 
for the process being considered. Thes amplitudes are c losely 
related to the functions p,„(x), and qm{x) in equation ( 1.19 ). 
The effective potential l/°*(x) experienced by the vapor as a 
result of the existence of a condensate and its quasiparticles. 



V^i^c) = 2wp(x) 



(1.25) 



2MiV|^jv(x)|^ 



{n™|p™(x)|2 + (n„ + l)|q„(x)|2} (1.26) 



where the mean occupation number of th e quasiparticle rn is 
defined in terms of the operators of (1.19) and the condensate 
band density operator pc- 

fim = Tr {bl^bmPc) , (1-27) 
and is time-dependent if pc is time-dependent. 

2. Condensate-band equations of motion 

The methodology used does not demand that any further ap- 
proximation be ma de to t he condensate-band effective Hamil- 
tonian as given in (1.13). However the interactions with the 
non-condensate band are most natur ally ex pressed in terms of 
quasiparticles, as explained in Sect. [ B 3 , and this will only 
make sense when the condensate-band Hamiltonian is ex- 
pressed in terms of the Bogoliubov theory. To a first approxi- 
mation it is reasonable to treat the condensate-band as a con- 
densate plus a gas of non-interacting quasiparticles, and this 
is all we do explicitly in this paper. 



There is no difficulty in principle in including the interac- 
tions between the quasiparticles which would arise in a full ex- 
pansion of the condensate-band Hamiltonian using the quasi- 
particle states as a basis. However a practical method of do- 
ing this would probably depend on the specific problem under 
consideration. 



II. USE OF WAVELETS FOR FOR TRAPPED ATOMS 

In QKI we introduced a wavelet expansion of the field op- 
erators. The wavelets are a complete set of orthonormal one- 
particle wavefunctions which are only significantly different 
from zero in a phase space volume around a point TiK, r 
in phase space. By this we mean that the wavelet function 
vk(x, r) is only significantly different from zero in a co- 
ordinate space volume tt^ / A"* around r, and that their expres- 
sion in terms of momentum eigenfunctions with eigenvalue p 
is only significantly different from zero in a momentum space 
volume (2?;,A)'^ around p — ?iK. This expansion enables one 
to express the quantum states of the Bose gas in terms of the 
numbers of particles associated with the wavelet K, r; which 
we can loosely think of as the number of particles in a phase 
space "cell" of volume at the point K, r. This gives a fully 
quantum mechanical description of the physics, but also con- 
nects very directly the classical description of the particle as 
having position r, and momentum KK. Using the descrip- 
tion, we can talk about scattering between different momen- 
tum values (the Boltzmann collision term), and flow through 
phase space (the Boltzmann streaming term). Because it is 
fully quantum mechanical, it is reasonably straightforward to 
use this phase space description for the non-condensate band 
i?Nc at the same time as we use a description in terms of en- 
ergy eigenfunctions in the condensate band Rq. 



A. Formulation of wavelets 

1. Wavelet properties 

The details of the construction of the wavelets are contained 
in Appendix^. The essential properties of the wavelets which 
we shall use here are as follows. We write the one-particle 
trapping Hamiltonian as 



'one particle) 



p 

2to 



y(x) 



(2.1) 



we can choose a set of wavelet functions (x, r) with the 
following properties. 

1 . The wavelets are a complete orthonormal set 
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J (i^xwQ(x,r)'u;Q/(x,r') = (5QQ'(Srr', 
^w*Q{x,r)wQ{yi',r) ^ (5(x-x'). 



(2.2) 
(2.3) 



2. The wavelet function ?i'Q(x, r) is only significantly differ- 
ent from zero when x w r. 

3. The mean energy of the wavelet is independent of r, and is 



E{Q) 



2m 



(2.4) 



The mean momentum associated with a wavelet is a vector 
?iK(Q, r) which is such that 



2m 



2m 



V{r) 



(2.5) 



The direction of Q is determined by the condition that the Li- 
ouvillian flow in phase space does not change Q; this means 
that K(Q,r) can be inverted to give a function Q(K,r), 
which can be regarded as a smooth function of continuous 
variable on a sufficiently large scale, and that this function 
satisfies the equation 



hK ■ Vr vy(r) • Vk 



Q(K,r) = 0. 



(2.6) 



Clearly, this is compatible with (2.5). We can fix the function 
by requiring that 



Vir) = 



Q(K,r)=K. 



(2.7) 



These requirements can be all satisfied simultaneously, as is 
demonstrated in Appendix 0. 



2. Field operator in terms of wavelets 



1 . The field operator can be expressed as 

i^(x) =^WQ(x,r)^Q,r (2.8) 

Q,r 

= ^V'q(x), (2.9) 
Q 



and 



(2.10) 



2. For r, Q sufficiently large that they may be regarded as be- 
ing almost continuous variables, the commutator of the V'q (x) 
functions has the asymptotic form 

[VQ(x),^t^,(x')] = <5q,q, ^u;Q,r(x)u;2j,r(x') (2.11) 



ffQ(x,x'). 



(2.12) 



where the function gQ(x, x') is only significantly different 
from zero when x w x', and is normalized to 



j d3x'5Q(x,x') - j d^xgqlx.x') = 1. (2. 



13) 



The region of validity of this asymptotic form coincides es- 
sentially with that of the WKB approximation for the wave- 
functions. 



3. Expression of the many particle Hamiltonian 



1 . The many particle trapping Hamiltonian — that is, exclud- 
ing interactions (13) — can be expressed as 



r Q 



EE^(Q''-''-')^Q.r^Q.r'. (2.14) 
rr' Q 



2. The object A^(Q, r, r') represents the quantized version of 
the Liouvillian motion, and defines a quantity A/ (K, K', r, r') 
through 

A/(K(Q,r),K(Q,r'),r,r') =A^(Q,r,r'), (2.15) 



a definition which is only valid for (r, K), (r', K') on the 
same Q surface. When applied to a sufficiently smooth func- 
tion i^(K,r), the M function generates the Liouvillian mo- 
tion through 



Y, M(K,K',t,t')F(K',r') 

fMC^_W(^| (2.6) 

y m n ] 



K 



3. T he interaction part of the many-particle Hamiltonian, Hi, 
(1.3) can be written in the form 



Q2,Q3,Q4 



(2.17) 



1234 



in which 



^Qi,Q2,Q3,Q4 

= Jd'^l rf=^x'^.(x-x')V'^j,(x)<(x')^Q3(x')V'Q4(x) 

(2.18) 
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B. Definition of the condensate band and non condensate bands 

The formulation in terms of wavelets so far has been ap- 
propriate for a situation in which no condensation occurs. We 
now consider the situation in which there is significant con- 
densation. There are two aspects to be considered: 

1. We must define the condensate-band and the non- 
condensate-band in terms of the Q-band formulation. The 
condensate band, which we will call Rc, will consist of all 
Q with energy less than a value _Er, that is all Q such that 
?i^Q^/2m < Sr. The non-condensate band, which we will 
call i?NC. consists of the complement of Rc', that is, all states 
for which ri^QV277i > Er. The non-condensate-band will 
be described in terms of wavelets, while the condensate-band 
will be described in terms of many-body eigenfunctions. 

2. If there is significant condensation, there will be significant 
mean-field effects, which will give rise to a cor rection to the 
the trapping potential, as discussed in Sect. [B 2. For the non- 



condensate-band this means that the definition of the wavelets 
used will be changed, because these are defined in terms of a 
particular potential, as in (2.1). 



1. Commutation relations 

The procedure we will now follow is to resolve the field 
operators into Q-bands only within the non-condensate band, 
where we will find that the Q-bands are effectively thermal- 
ized, and are thus uncorrected with each other In the con- 
densate band, significant correlation effects may arise because 
the interactions with the condensate are not negligible in this 
band — this leads to the use of a single field operator the de- 
scribe the behavior of the whole of Rq ■ 

We now define the non-condensate band field operator 
V'Nc(x), the field operator restricted to i?NC> based on the 
Q-band form for the field operator ( 2.8 , 0), as 



NC 



V'Nc(x) = ^^«;Q(x,r)AQ.r (2.19) 

Q r 

NC 

Q 



(2.20) 



where the notation X^q*^ means the summation over Q e 
i?Nc> with a corresponding notation for summations in Rc- 

The condensate band field operator will be called (/)(x), and 
is defined by 

0(x) = 7/.(x)-Vnc(x) (2.21) 
The commutation relations for these operators are 

Q r 

= <?c(x,x'). (2.22) 



The (f> and tpQ operators commute with each other, an d the 
commutation relations of the iPq{x) are as in ( 2.1 1 , 2.12 ) 



2. Separation of condensate and non-condensate parts of the full 
Hamiltonian 



We can now write the full Hamiltonian in a form which sep- 
arates the three components; namely, those which act within 
i?NC only, those which act within Rc only, and those which 
cause transfers of energy or population between Rc and i?NC ■ 
Thus we write the trapping plus kinetic parts of the Hamilto- 
nian as the sum of three terms 



-ffkin + Ht — Ha + Ht, + Hc,l 



(2.23) 



in which the part which gives the average energy of the non- 
condensate atoms in each Q-band is 



r Q 



(2.24) 



The part which represents Liouvillian transport of the non- 
condensate atoms in each Q-band is 

NC 

i?6 = EE-^(Q'^'^')^Q,r^Q>r'. (2.25) 
rr' Q 

The part which represents the non-interaction part of the 
condensate-band Hamiltonian is 



HcA = Jd^xcP^K) 



2m 



+ yT(x)U(x). (2.26) 



The interaction part Hi as defined in ( |1.3| ) can now be resolved 
into: 

1 . The part involving only (f) operators, which is the self inter- 
action within Rc, which we call 



H 



(2.27) 



2. The part involving no cf) operators, which gives rise to scat- 
tering of the particles within i?NCj which we call 

Hl^NC = 

d^x J rf^x- ""^^^ ^'^ <cW<c(xO^Nc(x)V'Nc(xO. 

(2.28) 

3. There are finally the terms involving operators from both 
bands, which cause transfer of energy and/or particles be- 
tween Rc and i?NC- We call the parts involving one op- 
erator 
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(1) 
i,c 



d^^ / d3x'u(x-x')<cWV^Nc(x')'/'Nc(x)<^(x') 



+ j d^^j d'x'w(x-x')0t(x)V'J,c(x')V'Nc(x)VNc(x'). 

(2.29) 

The parts involving two operators are called 

rr(2) _ 
-"l,C = 

d^x / d-''x'M(x - xO^'Jjc (x)-/-^ (x')^NC (x)</)(x') 



+ / d^x / d3x'^(x-x')</'^(x)Vj,c(x')V'Nc(x)0(x') 



^7^l,cWV'ic(x')'/'(x)0(x') 



d3x / rf3^,!4^i_iO0t(x)^t(x')^j,c(x)V'Nc(x'). 



(2.30) 



The parts involving three (f) operators are called 



(3) 
I,C 



d^x / dVu(x-x')<?it(x)<?it(x')0(x)V'Nc(x') 



+ I d^yi (i3x'y(x-x')^/'J,c(x)0^(x')'/'(x)0(x'). 



(2.31) 

Our notation is not entirely systematic, since -ffi.NC involves 
only the non-condensate-band operators, whereas iJi.c in- 
volves operators from both the condensate-band and the non- 
condensate-band. 



3. Mean field corrections 

In practice the thermalized part of the condensate-vapor 
system consists of nearly all of the levels above the conden- 
sate level, and it is advantageous to include this fact explicitly 
as early as possible in the analysis. Thus, we would like to 
put the division between the condensate and non-condensate 
bands at as low an energy as possible. The main difference 
between the condensate band and the non-condensate band is 
the nature of the excitation spectrum. It is relatively simple to 
include corrections to the energy levels of the non-condensate 
band provided these energy levels are particle-like — that is, 
they are states with a definite energy and a definite number of 
particles. In fact, computations have shown only a small pro- 
portion of the excitation spectrum is not-particle like, and this 
is at quite low energy [|5]|. However, the energies of lower 
energy particle-like excitations are affected by the presence of 
the condensate and of each other. We therefore must consider 
mean field effects. There are three kinds of mean field effect 
which can be explicitly included: 



1. The average effect of the condensate band on the non- 
condensate band; 

2. The average effect of the non-condensate band on itself; 

3. The average effect of the the non-condensate band on the 
condensate band. 

Our aim will be to estimate these average effects explicitly, 
and include them in the two Hamiltonians used to describe the 
two bands. The sum of these two terms is then subtracted from 
the interaction part of the Hamiltonian, which is preserved as 
a correction term. The description is thus in principle exact, 
since the total Hamiltonian is not changed — only the basis for 
the perturbation theory is changed. 

If we suppose that p{x) is a c-number estimate of the 
condensate band density, and n(x) is a c-number estimate 
of the non-condensate band density, then the mean field ef- 
fects 1. and 2. give an equivalent potential acting on the non- 
condensate 



T/mf 



(x) 2u(^p{x) +n(x) 



(2.32) 



while the effect 3. is given by an equivalent potential acting 
on the condensate 



(2.33) 



The mean field effect is included by adding and subtracting 
a term 

J d3x{FN"g(x)4cWV'Nc(x) + Fc"^f(x)0t(x)0(x)} . 

(2.34) 

Effectively, this changes the potentials used in defining the 
two bands. For the non-condensate band, the change is 



(2.35) 



with a consequent redefinition of wavelets. Thus the 
non-conde nsate band Hamiltonian retains the form as in 
( 2.24 , 2.25 ), but the wavelets are now defined in terms of the 
effective potential V^q{x.). 

For the condensate the change is 



which now transforms the He to H^^, defined by 



(2.36) 



2m 



+ yc'*f(x)^ 0(x)+i?c,: 



(2.37) 



The extra terms are then subtracted off the interaction part of 
the Hamiltonian, giving 



(2) 
I,C 



(2,eff) 
I,C 

(2) 
I,C 



d3xV^N"c(x)V'J,c(^)^Nc(x) 



d3xV(i"f(x)?i(x)0t(x). 



(2.38) 
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We use the term fff ^ for the sum of all of these interaction 
terms; 



rrcff 



(1) 
i,c 



(2,off) 
I.C 



H 



(3) 
LC- 



(2.39) 



The description is still exact, since the mean-field terms have 
been added and subtracted, but an estimate of the mean effect 
of the interaction has been explicitly included in the basis de- 
scription of the condensate and non-condensate. T he method - 
ology does not depend on the particular choices of ( 2.32 , 2.33 ). 
It is possible, for example 

1. To choose the mean field potentials to be fixed quantities; 

2. To choose them to depend on the numbers of particles in 
the condensate and non-condensate bands; 

3. To include an explicit time dependence. However, if this 
is done, the derivation and form of the master equation in 
Sect.|ll| would be altered. 

The choice of the estimates for fi(x) and p(x) is not made 
at this stage, but will be done as part of the derivation of the 
master equation, since their choice must be such as to opti- 
mize the basis in which the master equation is described. 



III. DERIVATION OF THE MASTER EQUATION 

Using the separation into condensate band and non- 
condensate band operators we can now proceed to develop a 
master equation in a a way that combines the methodologies 
ofQKland QKIII. 



A. Projectors in the Q-bands 

We can now introduce the definitions of the projectors 
which we shall need, in much the same way as in QKl. The 
states of the system will be written in the form 



*c,n) 



(3.1) 



in which represents the condensate degrees of freedom, 
and n is a vector of integers nq.r representing the numbers of 
atoms in the Q-cell Q, r. We will want to consider a projector 
definition of the form 

PN|*c,n) = |*c,n) if ^n(Q',r) = 7V(Q), 

r 

= otherwise. (3.2) 

This projector projects only in i?Nc- leaving the condensate 
state description 'I'c completely unaffected. 

The projectors on the density operator are defined as in QKl 



= wn. 

The complementary projector Q is defined by 



(3.3) 
(3.4) 



Q = l 



(3.5) 



N 



where, by the definition ( [3.2| ) of V-yi, the range of the index N 
includes configurations entirely within i?Nc . 

Thus the terms v-^ will be the parts of the density operator 
for which we can derive a master equation, while the comple- 
mentary part 



w 



Qp 



(3.6) 



is the part which is neglected. 

The quantities are thus density operators in which N 
represents the Q-band configuration of non-condensate-band 
particles, and in which no restriction of any kind is placed on 
the configuration of the particles in the condensate band. 



B. Formal derivation of the master equation 

The formal derivation of the master equation follows much 
the same methodology as in QKl and QKIII. In the present sit- 
uation, for the purposes of developing the master equation we 



can divide the operators into the various parts as in Sec. II B 2 



but also taking account of the fact that we use the effective 
Hamiltonians which include estimates of the mean field ef- 
fects: 



rrcff zjcff 

- ^C,l ' 

rrcff crcff 

-"I — ^I,NC 

rrcff rrcff 



C,2 



HI 



i,c 
s 



[see (2.37)], 
[see ( pl , ^39| )] 



(3.7a) 
(3.7b) 
(3.7c) 



In the last of these, we understand H^^ and H?^ to be con- 



structed in the same way as ( |2.24| , [2.25| ), b ut us ing the effec- 
tive potential V^q{x) innstead of ^(x) in (2.1) as a basis for 
counstructing the wavelets. 
There are the relations 



[Hf,V^p]=0 



(3.8a) 
(3.8b) 
(3.8c) 



which follow directly from the construction of the operators 
and projectors. 

The equation of motion is 



-[HI 



cff 



Hf + Hf 



{Ca+Cb + Cl+Cc)p. 



(3.9) 
(3.10) 



Notice that, for compactness, we do not explicitly write C^^ 
for the Liouvillians. 

From this form it is straightforward to write the master 
equation in the same form as in QKL As in QKl, we use the 
Laplace transform notation for any function J{t) 



(3.11) 
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Using this notation and the relations (3.8a-3.8c), the master 
equation takes the form 



SUN (s) - VN (0) = £h€'N (s) + "Pn-C/ ^ WM (s) 

M 



M 



(3.12) 



In this form the master equation is basically exact. We shall 
make the approximation that the kernel of the second part, the 
[ term, can be approximated by keeping only the terms 
which describe the effective Hamiltonians within Rq or i?Nc- 
namely the terms Ca and Cc- Thus the master equation we 
shall use is written 



sun(s) - wn(0) = Awn(s) + V-hCi ^ vm{s) 

M 

+V^Ci[s- Ca- Ccr^QCiY,im{s) (3.13) 

M 

This is essentially the same approximation as in QKI, apart 
from the inclusion of the condensate band self-interaction 
term in Cc- However it should be noted that Rq is a rather 
wide band compared to the the K = band of QKI, and 
contains significant dynamics, which we do not assume to be 
thermalized — rather, we expect these levels to become ther- 
malized as a consequence of the equations of motion. 



d3xyiJc'(x)V'Ij(x)VQ(x) 



(3.14) 



In (3.14) one can see clearly the "Hartree" term in the third 



line, and the "Fock" term in the fourth line. However in the 
remainder of this paper we shall use the delta function form 
of the interaction potential ■u(x — x') uS{x — x'), and the 
difference between these terms and the second and third lines 



of (2.30) disappears. None of our results will depend heav- 



ily on the validity of this approximation, but the formulae do 
become considerably more compact when it is made. 

We thus arrive at the final form of the forward scattering 
term 



off 
'forward 



Qi#Q2 Q I 



QQQ 



2u / d3xV'^j(x)VQ(x) {0t(x)0(x) - p(x) - n(x)} 



-2u / d3x,l(x)(/|t(x)0(x) 



(3.15) 



D. CoUisional terms 



C. Reversible terms 

1. Streaming Term 

This arises from the term Cj^v^ ( t), an d takes an essentially 
classical form, based on of ( 2.14 - 2.16 ), in the degree of ap- 
proximation we shall later use. 



2. Mean-field and forward scattering terms 



These arise from the term Pn-C/ ^ wm(s). Using the 

quantity Z^Qi,Q2,Q3,Q4 as defined in ( 2.18| ), we find that these 
terms are equivalent to a Hamiltonian term of the form 



The remaining parts of the master equation (3.12) lead to 



Q2,Qi,Q2 +^Qi.Q2,Q2,Qi 



Qi#Q2 
Q I 

+ Jd^^J dVu(x-x')Vj3(x)0t(x')</)(x')?AQ(x) 
+ [ d^x /"dVii(x-x')V'ij(x)<^t(x')0(x)^,Q(x') 



genuinely irreversible collision terms, both within the non- 
condensate band, and between the two bands. 



1. Collisions within the non-condensate band 

The interactions which involve only particles in i?NC can 
cause collisions between the particles in the noncondensate 
band. Thus we consider the parts of the second line of ( 3.13 ) 
which arise from the part 



^nA,nc[s -Ca- Cc]-^QCi^Nc ^M^'*)- (^.16) 



M 



The result of this process can be seen by considering a par- 
ticular term in the expansion of (3.16). The inversion of the 



Laplace transform would give a term like 



/>oo 

U{12M) / dr cxp{{Ca + Cc)T}U{4'3'2'l')vN{t - r) 
Jo 



(3.17) 



The only surviving term after projection will have 

1 = l',2 = 2' or 1 = 2',2 = 1' 
and 

3 = 3',4 = 4' or 3 = 4', 4 = 3' 
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and we can write this as 

Z^(1234) / dTe-'^'^^-''*''U{A321) 
Jo 

where, as in QKI, 

AcJi234 = LOi + UJ3 — UJ2 — UJl, 

but where we use the notations 



hu), 



(3.18) 



(3.19) 



(3.20) 
(3.21) 



We now make the usual approximation that, in the time inter- 
val over which the r integration is significant, we may ignore 
the irreversible part of the time evolution of wn, so that the 
exponential terms simply set t — r — ^ t, and we get simply 



/>oo 

Z^(1234) / dre" 
Jo 



-iAwi234T 



U{A321)vT<s{t). 



(3.22) 



which leads to the terms in the QKME which describe the 
scattering of the modes in i?Nc in exactly the same form as in 
(QKl;75d,e). 



2. Interaction between Rc and Rnc 

The treatment of these terms depends on writing the in- 
teraction Hamiltonian terms in a form which uses wavelet 
functions for the non-condensate-band operators, and energy 
eigenfunctions for the condensate band trap operators as fol- 
lows: 



rr(l) 



(2, off) 
I,C 



(3) 
I,C 



^ / d3xZ432(x)^Xt„(x)+h.C. 

234 m 

2^ / rf^xZ42(x)^C/t(x)+h.c. 

42 r 

f rf^xZ4(x)^v;t(x) + h.C. 



(3.23) 
(3.24) 
(3.25) 



The non-condensate band is represented by the Z operators; 
an abbreviated notation for 



^432 (x) = U?/'4(x)V'3(x)-0|(x), 
Z42(x) = li?/'4(x)'02(x), 
Z4(x) = U?/'4(x). 



(3.26) 
(3.27) 
(3.28) 



The condensate band operators are defined by 

</)(x) =^X,„(x) (3.29) 

m 

0(x)0t(x) =^C/,(x), (3.30) 

r 

0t(x),^(x)0(x) =^K-(x). (3.31) 

r 

The X, U, V operators are, as in QKIII (where they are dis- 
cussed more extensively), eigenoperators of the condensate- 
band effective Hamiltonian, with eigenvalues defined by 



[Hf,Xra{^)]=-hemXrn{^)., 
[Hf,Ur{^)]^-hl3rUr{^)., 

[if£«,F,(x)] =-?i/3;K(x) 



(3.32) 
(3.33) 
(3.34) 



In the Bogoliubov approximation, these operators have rel- 
atively simple expressions in terms of quasiparticle creation 
and destruction operators — the precise relationship is given 
in Sec.llD of QKlll; especially (94-99) which gives the ba- 
sic definition, and Sec.llD3 which gives the expression of the 
more complicated Ur{x) and Vr{x.) operators. 

We omit terms involving 0(x)(/)(x), or its Hermitian conju- 
gate, since these are always very non-resonant. We use the ro- 
taing wave, or random phase approximation as in QKlll, and 
derive a master equation which is a blend of the formalism of 
QKI and that of QKlll. 



E. Master equation terms 

Putting these all together we get the contribution to the mas- 
ter equation which describes coupling between Rq and i?NC 
so that the full quantum kinetic master equation now takes the 
form 

i'N(i) = i'N(i)|c + WN(i)|NC + WN(i)|MF + i'N(i) I coupling 

(3.35) 

in which the various terms are defined as follows. 



1. Condensate band term 



This is 



(3.36) 



where, as in (2.37) 



= J d^x (/)t(x) f^- 



^ + y^glx)) 0(x) + ^ / d3x0t(x)<^t(x)0(x)0(x) 



(3.37) 
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2. Non-condensate band term 



■i'N(i)|NC = ~^ [i?6,WN(i)] 



^ [i?F,NC,fN(0] 



+ ^^i5(Aw(e)){2Z^(e)wN-e(i)Z^^(e) -Z^^(e)W(e)uN(t) - WN(t)Z^t(e)Zi(e)}. 

e 

Here we define the quantity 



(3.38a) 
(3.38b) 
(3.38c) 

(3.38d) 
(3.39) 



Q 



Qi7^Q2 



Q 



which is that part of the forward scattering term ( 3.14 ) — including the mean field correction — and affects only the non- 
condensate band. 

As in QKl we have defined the vector in the space of M, N, called 61234 which can be written 



61234 


= 61234 


— —61234 


— —61234 


= 1 




Qi 


Q2 


Q3 


Q4 



(3.40) 

with all other components are zero. Thus 61234 represents a collision, and we have used the notation as a shorthand for a 
summation over Qi, Q2, Q3, Q4. 



3. Mean-field coupling term 



This is the term which gives rise to mean-field effects, and arises from that part of the forward scattering term (3.14) which 
mixes together the condensate-band and non-condensate band operators. It is 



with 



ff^f = 2uj d^x I (^^^.tj(x)V^Q(x) - n(x)^ (^0^(x)0(x) - p(x)) 



(3.41) 



(3.42) 



This term is a relatively small coupling term, since the major mean field effects have been canceled. For ease of writing, a c- 
number term 2 J (Px n(x)p(x) has been added, leading to the factorized form given. Of course this does not affect the equations 
of motion. 



4. Coupling terms 

These terms give all transfer of energy and population between Rc and Rj<sc, and are: 



coupling 



^ y d^X y d^x' ^2gX]{'^''(^'^234-ern)^432(x)Zt32(x')Xt„(x)X,„(x')«N(i) 



+6P{- AUJ234 + e™)«N(i)^432(x)Zl32(x')Xt„(x)X™(x') - 25(AtJ234 - e™)X™(x')^l32(x')«N-e,32(i)^432(x)Xt„(x) 

(3.43a) 

('5''(-^^234 + 6,„)Zi32(x)Z432(x')^r„(x)Xi(x')«N(i) 



432 rn 
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(t)Zl32(x)Z432(x')X™(x)Xt„(x') - 2(5(At^234 - e™)X (x')Z432 (x')z;N+e,3. (i) ^l32 (x) 

(3.43b) 

Z42 (x) (x' ) f/^ (x) [/r (x' ) wn (i) 

42 r 

+5P(-At^24 + /3.)wN(t)^42(x)Zl2(x' )C/t(x)[/r(x') - 2(5(Aa'24 - /3r)f^r(x')^l2(x')'^N-e42 (i)Z42(x);7t(x)| (3.43c) 

+ |jP(^^ - /?;)^4(x)zi(x')F;(x)K(x')t'NW 

4 ,. L 

+SP{^uji + /3;)z;N(i)^4(x)Z](x')F;(x)l/,(x') - 2(5(^4 - /3;)K(x')^l(x')«N-e4 W^4(x)F;(x)| (3.43d) 

+ E E ('^''(-'^4 + /3;)zi(x)Z4(x')K(x)F;(x')m(i) 



4 r 



+SP{oj4 - /3;)t;N(i)^4(x)^4(x')F.(x)F;(x') - 2(5(cj4 - /3;)K^(x')^4(x')«N+e, (i)^l(x)V^.(x) 



(3.43e) 



These expressions use the function 

SP{uj) = 6iuj) - --. 



(3.44) 



However, for simplicity we shall usually neglect the dispersive 
effects which arise from the principal valu integral by setting 
^ 6 everywhere. 



5. Full master equation 

The full density operator, in the approximation being used 
for these equations of motion, can be written as 



(3.45) 



N 



where the summation is over all N in the non-con densate 
band. Summing over N in (3.36), (3.38a-3.38d),and in 



( |3.43a| -3.43e), it is easy to see that one obtains an equation 
for p(t) which is of exactly the same form, but in which 
■yN(i),WN-e(i) p{t) forany N,e. 



IV. INTERPRETATION OF THE MASTER EQUATION 

The master equation gives a representation of transitions 
due to collisions between states, defined by wavelets in the 
non-condensate-band, or eigenfunctions in the condensate 
band in their respective effective potentials which already in- 
clude the majority of mean-field effects. The collisional terms 
include those purely within the non-condensate-band as in 
(3.38d), and those which give rise to energy and matter trans- 
fer between the two bands, as in (3.43b-3.43e). In addition 



to the truly irreversible collisional terms, there are also terms 
which arise from principal value integrals, which cause energy 



level shifts, since the y take the form of a commutator. These 
terms are explicit in (3.38c), and are implicit in the terms 
in (3.43t-3.43e), and are analogous to the Lamb and Stark 



shift terms which arise similarly in quantum optics. They are 
equivalent to energy level shifts found by using second order 
perturbation theory. 



A. Stationary solutions of the master equation 

The residual mean-field terms and the level shifts are not 
expected to be large, nor do they in any way change the es- 
sential picture of transitions. To the extent that we can neglect 
these terms, the quantum kinetic master equation in the form 
(3.35) conserves 



off 



rrcff — Efoff 
-"c = J^tot 



A^NC + iVc = ^, 



(4.1) 
(4.2) 



so that any function of these provides a stationary solution. 
The grand canonical stationary solution takes the form 



p cx cxp 



H. 



cff 



kT 



) exp 



kT 



(4.3) 



but the fact that A^nc + Nc is conserved means that the re- 
striction of the grand canonical form to any fixed value of 
A^NC + A'c must also be a solution. The correlation functions 
of the non-condensate operators are quite straightforward 



r.r'OQ.Q' 



7Vt(Q) 



(4.4) 



with 
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Nt{Q) 



exp 



1 -1 



exp 



kT 



- 1 



(4.5) 

-1 

(4.6) 



We can also deduce that the stationary field correlation func- 
tions can be written 



(vij(x)^Q,(x')) 



Sq.Q'Nt{Q) WK(Q,r),r(x)wK(Q,r).r(x') (4.7) 



;k Q, 



• (x-x'; 



5q(x,x'), (4.8) 



where the last line comes from the asymptotic form (2.12|. 
(These forms will also be valid to a very good degree of 
approximation even in a canonical ensemble formulation in 
which A^NC is fixed.) 



and when we are dealing with the stationary state, this is given 
through (0) by as 



n(x)= J2 ^T(Q)ffQ(x,x) 

QS-Rnc 



£(K,x)>_Eh 



Here 



i?(K,x) = ^ + yi5g(x) 



2 m 
2m 



(4.10) 



(4.11) 



(4.12) 



+ 1/t(x) + 2u(p(x)+7i(x)). (4.13) 



This amounts to a local density approximation for the atom 
density of the non-condensate-band. 



B. Relationship to tlie Hartree-Fock-Bogoliubov metliod 



2. Condensate-hand density p(x) 



The inclusion of the mean-field terms defines the stationary 
solutions of the Hamiltonian part in a form which is closely re- 
lated to the Hartree-Fock-Bogoliubov method |20|. The task 



which must be accomplished here is the construction of best 
estimates of the condensate-band density p(x) an d of the non- 
condensate-band density n(x) as defined in Sect. [IB 3. 



1. Non-condensate-band density n(x) 

The non-condensate-band density must be chosen in the 
form whi ch tak es account of the projected form of the inter- 
action in (3.14), i.e.. 



n(x)= (^Ij(x)^q(x)). 



(4.9) 



It is now necessary to specify p(x), which means we need 
to know the condensate state. From {\3) we see that this 



requires t he di agonalization of iJg* — ^Nq, where H^^ is 
given by (|37j). This can be done approximately by using the 
Bogoliubov method [^, generalized to its number conserving 
form as given in [0], but using the V(^*(x) potential, which in- 
cludes the mean field terms arising from the non-condensate- 
band. This would mean that we write 



(x)«AUiv^(x 



1 



=x(x) , 



(4.14) 



where the operator A is a destruction operator for the number 
of atoms in the condensate-band, so A'^ A has the eigenvalue 
Nq, and ^TVc (x) is a solution of the Gross-Pi taevskii equation 
for Nq particles in the effective potential. Thus 



2y72 



2m 



<Nc{^) + K?''(x)eJVc(x) + uNc |eiVc(x)|'eA'c(x) = A^Wc^JVclx). 



The residual part x(x) is determined in terms of quasiparticles, thus 

X(x) = ^ [p™(x)fe„, + (7,„(x)6j„ 

where Pm(?^) ^nd gm(x) diagonalize the residual Hamiltonian 

= d?^x\^)^\{^) + I rf^xxnx)^c'*'(x)x(x) 



uNc 



{^Na{^)xH^)Y + ^(Gc Wx(x))' + xt(x)x(x)(2u7Vc|^^c(x)|' - /i^c) 



(4.15) 



(4.16) 



(4.17) 
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The diagonalization of 7^3 is equivalent to solving the 
Bogoliubov-de Gennes equations with a potential 

Vq^(x), and with an upper cutoff energy equal to Eji — /xjv^.. 
The wavefunctions p,„ (x)and (x) are the projections of the 
resulting solutions onto the subspace orthogonal to ^jvc (x)- 

This procedure gives solutions with any number Nq, ^Vnc 
of condensate- or non-condensate-band particles. The value 
of the /?(x) which arises is then 

p(x)«7VcieA'c(x)p 

+ ^ {"-mbm(x)|" 



(n™ + l)|g™(x)|2}. (4.18) 



elude the master equation coupling terms that the condensate- 
band and non-condensate-band are forced to evolve to an equi- 
librium in which the chemical potentials and temperatures are 
the same. The relevant equations of motion can then be com- 
puted in various degrees of simplification, and it is these that 
should be compared with the time-dependent Hartree-Fock- 
Bogoliubov equations of motion. 

5. The energy levels, wavefunctions and densities given by 
this process are to be seen as a way of choosing a good ba- 
sis in ter ms of which the processes described by the master 
equation (3.35) can be described. As well as the irreversible 
processes, there will be corrections to these which also arise 
from the master equation, and which can be in principle com- 
puted. 



C. Validity of tlie metliod 



1. The Bogoliubov approximation used here will be valid 
when the majority of the atoms in the condensate-band are 
in fact in the condensate level. This is not normally the ma- 
jority of the atoms in the system, since there may be a very 
large number of atoms in the non-condensate-band, which are 
treated essentially by a local density approximation. 

2. The values of n(x) and p(x) are used then to determine 
the mean field potentials, and hence the appropriate energy 
levels which provide a basis within which the master equation 
now describes time evolution, including transfer of energy and 
atoms from one band to the other, as well as the residual re- 
versible processes. 

3. The procedure for determining this basis is essentially 
equivalent to the Hartree-Fock-Bogoliubov method in the 
Popov approximation, but only in its time independent form. 

4. Note that the chemical potentials and temperatures of the 
stationary solutions of the construction so far used can be de- 
termined independently of each other — it is only when we in- 



V. TIME-INDEPENDENT LOCAL EQUILIBRIUM 
APPROXIMATION 



In this section we will show how the methodology of QKIII 
arises. We assume that the scattering in _Rnc is strong, and 
we are only interested in the behavior of the condensate and 
the modes in Rq. We then factorize the total density oper- 
ator into a quasi-thermalized non-condensate part ps, and a 
density operator pc{t)^ which describes i?c- Thus 



pit) = PB® Pc{t), 



(5.1) 



N 



and is of course time-independent. We use a local equi- 
librium approximation, in which ps is such that the averages 
of the products of field operators in i?NC are quasi-thermal, 
in much the same way as was used in the derivation of the 
Ueh ling-Uhlenbeck equation in QKI. Thus we write, as in 



(5.2) 



This allows for thermalization locally, in contrast to (4.8), for which Fq ^ ^"^^ ^ Nt{Q) which give a global thermalization. 

We use the Gaussian factorization properties of the the averages of the products of more than two field operators, so that for 
example, we can write 



(V'qi(x)?/'q2(x)?/'q'^(x')V'q;,(x')) = (^Qi,q'/q2.q;, + Sq^.q'Jq^^'J 9q{^ - y^'fFQ, ^ 

^^{,K(Qi,2£±^)+iK(Q2,^)}(x-x')_ 



(5.3) 



A. Master equation 



Tracing over the non-condensate modes, we get 



Pcit) 



0t(x) 



(x) , PC 



(5.4a) 



16 



+ I S^j d^x' ( G'^Hx - X', £,„) (2X„,(x)pc^,Ux') - pcXl{^)X,-a{^) - Xl{^)Xm{^)pc) 

+ J2 G^^H^ - x', e™) (2Xt„(x)pc^™(x') - pc^™(x')Xt^(x) - X,„(x')Xt^(x)pc) 

m 

+ ^ A/(x - x', /?,) (2[/,(x)pcC/.t(x') - pcUU^')Uri^) - [/,^(x')C/r(x)pc) 

r 

+ E^+H^ - x', /?;.) (2V;(x)pcK^(x') - pcK^(x')K(x) - F;(x')K(x)pc) 

r 

+ Y E^~\^ - x',/?;) (2K.^(x)pcK(x') - pcV^.(x')y;(x) - Vr{^')Vj{^)pc) ^ 
in which the quantities E^^\ M, G^^-* are given by 

G'+^(x-x',w) = I^Trs i +Ci;2 -t^3 - w)Zj23(x')Zi23(x)pB i 

I 123 J 

G("^(x-x',w) = Ir^rs I Y^^^^ +u}2-i^3 - uj)Zi23{x)Zl2^{x')pB \ 

I 123 J 

M(x-x',cj) = ^Trs -u;2 -c^)^l2(x')^i2(xw| 

i?(-)(x - x',c.) = ^Trs |^<5(a;i - c.)Zi(x)Zt(x')pi3 1 . 



(5.4b) 
(5.4c) 
(5.4d) 
(5.4e) 

(5.4f) 



(5.5a) 
(5.5b) 
(5.5c) 
(5.5d) 
(5.5e) 



The master equation ( 5.4a - 5.4f ) is of the same form as given 
in QKIII, apart from 

1 . The expression of the non-condensate band summations in 
terms of Q-bands, rather than the continuum approximation 
chosen in QKIII as integrals over K variables. In this degree 
of approximation, this leads only to a change in the range of 
K as a function of x. 

2. The potential defining the Q-bands is the effective potential 
T^c(x), which depends on the occupation of the condensate- 
band, rather than simply the trapping potential as in QKIII. 



B. Interactions between quasiparticles 



As written in ( |337l ) and (| .4a[ ) the condensate-band Hamil- 
tonian is treated exactly. Using the Bogoliubov approxima- 
tion, we c an approximate the condensate-band Hamiltonian 
term (5.4a) in the form 



(5.6) 



Although the irreversible terms in (5.4a-5.4f ) give damping 
and thermalization, this approximation does omit the effect 
of the higher order terms, which can generate scattering of 



the quasiparticles by quasiparticles. Graham has pointed 
out that these can be significant, even though the quasiparticle 
density of states is rather low, essentially because the very low 
energy quasiparticles involve very large numbers of atoms. If 
this pertains, we have three possibilities. 

1) We can treat the quasiparticles in essentially the same man- 
ner as the particles in the non-condensate band; that is, we can 
assume they are thermalized, and compute scattering rates. 
There is no need, however, to use wavelets, since we can sim- 
ply use the discrete quasiparticle states, of which there should 
be relatively few compared to the non-condensate band states. 
This would correspond to Graham's methodology. 

2) It is possible that an exact treatment of the condensate band 
Hamiltonian could be feasible, provided it is realistic to take 
this as consisting of only a relatively few quasiparticle lev- 
els. Stringari [ |2^ ] has pointed out that in fact in the case of 
a trap, very few of the levels have any genuine phonon-like 
character. This could be done in combination with the use of 
a stochastic wavefunction method for the interactions with the 
non-condensate band. 

3) As an intermediate form of approximation, it could well 
be feasible to use a positive-P-representation methodology, 
rather like that used by Drummond and Corney [ |2^ ] for the 
condensate-band Hamiltonian in a quasiparticle basis. 
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We shall leave further developments of this aspect to an- 
other publication 



VI. TIME DEPENDENT LOCAL EQUILIBRIUM 
APPROXIMATION 

The local equilibrium approximation, when applied to the 
quantum-kinetic master equation of QKI was shown to give 
rise to the well known Uehling-Uhlenbeck equation [|l8|]. A 
similar method applied to this case leads to a description of 
the coupling of thermalized vapor to the condensate. Together 
with the equations of motion for the condensate density opera- 
tor as in Sec. 111.D.3 of QKlll, or the rate equation approxima- 
tions to these, this gives a description of the coupled motion 
of the condensate and non-condensate. From this description 
one can take further limits to get coupled equations of motion 
in the hydrodynamic limit, or in the weak-collision limit. 

To be precise, we define 



/k(x) = Fq(k,x)(x) 



(6.1) 



where F is as defined in (5.2). (This is the same as in QKI). 



Using the same procedures as in QKI, we can eventually ar- 
rive at the equation of motion: 



9/k(x) <9/k(x) 



dt 



9/k(x) 



9t uu 
5/k(x) 



dt 



+ 



dt 
9/k(x) 



mf 



dt 



9./k(x) 



dt 



(6.2) 



Of these, the first two terms represent flow and colUsions of 
atoms in the non-condensate-band, and the last three are irre- 
versible master equation terms, which represent the transfer of 
energy and atoms between the condensate-band and the non- 
condensate-band. 



The equation of motion (6.2) depends on the condensate- 
band state through the second term, which depends on 
the mean field determined by the mean condensate density 
p(x, t), as defined by (5.8). Thus we must complement (6.2) 
with an appropriate equation of motion for the condensate. At 
the simplest, we can simply use rate equations for the con- 
densate le vel o ccupationsand compute the condensate den- 
sity from ( 4.18 ). The condensate equations of motion should 
be wr itten in terms of the time-independent basis set as in 
( 4.16| ), which uses the time-independent mean-field potential 
Vq** (x). The deviation from the avera ge density n(x, t) in- 
duces the additional Hamiltonian term (3.41,3.42), which can 



also be expressed in the Bogoliubov basis. 

The process of condensate growth gives rise to an apparent 
time-dependence as follows. The basis set of states is differ- 
ent for every pair of values -/Vnc , — thus the set A^nc , 
can be regarded as defining a set of independent sectors, each 
with its own basis states. When the coupling terms between 
the condensate-band and the non-condensate band are taken 
into account, the effect is to transfer the system between sec- 
tors by the transitions A^nCj ~^ ^nc i 1, Nq T 1^ and 
to provide a new basis for the density operator. However this 
can be accounted by a slight modification to the Bogoliubov 
expansion — thus this apparent time dependence can in fact be 
fully accounted for by the existing procedures. 



A. Collisions and flow in the non-condensate 

1. Uehling-Uhlenbeck term 

This term is the part arising from streaming terms and col- 
lisions between atoms in the non-condensate-band with each 
other, and takes the form 



a/K(x) 



dt 



uu 



SK-Vx Vxl/i5g(x)- Vk 



m 



/k(x) 



-Rnc -Rnc -Rnc 



x{/k(x)/k.(x)[/k3(x) + 1][/k,(x) + 1] - [/k(x) + 1][/k.(x) + 1]/k3 (x)/k4 (x)} (6.3) 

Notice that the range of energies available is restricted to the condensate-band, which is consistent with the streaming term — thus 
there are no terms here which can transfer particles out of i?NC ■ 



2. Mean-field corrections to the Uehling-Uhlenbeck term 

The Uehling-Uhlenbeck term includes only the flow of 
the non-condensate which arises as a result of the (time- 
independent) effective potential V§q = 2up(x) + 2un(x) + 
Vt(x). However, when the condensate is not time- 
independent, the cancellation of the average values of Hmi as 
in (3.39,3-42) does not occur, since n(x) and p(x) are taken 



as time independent estimates of the mean values. Including 
these in the equation of motion for /k(x) gives rise to extra 
terms proportional to 



Q' 
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X (II v4'(y)'/'Q"(y) + 'P\y)^(.y) - ^(y) - p(y) 



(6.4) 



The commutator is non-zero only for Q' = Q since the com- 
mutators V'q(x) 6tc. are non-local. If we use the form (5.2) 
to evaluate this average this commutator does vanish — a non- 
vanishing result only appears when we allow corrections to 
this as a result of the non-equilibrium state implicit in the sit- 
uation. The detailed working is given in Appendix^ here we 
shall merely state the result, which is intuitively obvious, that 
for sufficiently smooth distributions in Q and x, we find 



a/K(x) 



mf 



dt 

= 2uVx {p(x, t) - p(x) + n(x, t) - fi(x)} • Vk/k(x) 

(6.5) 

where f7,(x, t) and p(x, t) are the time dependent versions of 
the condensate-band and non-condensate-band densities, i.e.. 



x, t) = 5] Tr U^, (x) V-Q' (x)ps (t)) (6.6) 



Q' 



p(x,t)=Tr((/.t(x)0(x)pcW). 



(6.7) 
(6.8) 



Thus the net effect on the UehUng-Uhlenbeck term is simply 
to make the replacement to the effective potential 



^Nc(x)-V^Ng(x,t) 

= Vt(x) + 2wp(x, t) + 2un(x, t), 



(6.9) 



that is, a time dependent version of the effective potential. 



3. Condensate growth term 

This is the principal term which gives rise to growth of the 
condensate. In the case that we use the Bogoliubov approxi- 
mation to the condensate, and the notation of Sec. lid of QKIII, 
to write this in terms of the two rate functions 



Rate^+^(K,x,w) = / d^Ka / d3K4 / d^k5(Atj2K4(x) - w)(5(K + K4 - Ka - k) 

(27r)2?i J J J 



Rnc Rnc 

X [l + /K(x)][l + /4(x)]/2(x)W,(x,k) 



(27r)2?i2 



Rnc Rnc 

x/K(x)/2(x)[l + /3(x)]W7(x,k) 



Rate[~^(K,x,cj) ^ , / d^Ka / d^K^ [ d^k6{Au;2K4i^) ~ uj)6{K + - K2 - k) 

(27r)2?i J J J 



Rnc Rnc 

x/K(x)/4(x)[l + /2(x)]W7(x,k) 



(27r)2?i2 



J f^^Kg y d^K^Jd^kS{Au;-K3i{x)^u;)S{K3 + Ki-K-k) 



Rnc Rnc 

X [1 + /k(x)][1 + /2(x)]/3(x)>V/(x,k) 



(6.10) 



(6.11) 



The term is then defined as 

9/(K,x) 



dt 



TVcRatcJ, '{K,ji,fiNc/^) - (^c + l)Ratc^^^(K,x,^Arc/?i) 
+ ^ J n„ Ratct,;) (K, x, i^lNc + e"c )/^) " + l)Ratc(+' (K, x, {^lNc + ^Nc)/^) 

m \ 

+ + l)Ratc^^^(K,x, {^inc - Etvc)/^) " rim Ratc^+^(K, x, (/^Wc - ejVc)/^) \ 



(6.12) 
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4. Scattering terms tlq = Nq (6.15) 

(6.16) 

We also need an abbreviated notation for the number oper- ^ — jy^n (6 17) 



ators of the condensate-band _ , , ^ 

(6.18) 

fii^iSlSi) (6.13) n^,n = Ncium + 1) (6.19) 

mi = {SiSl), (6.14) m_,„ (iVc + l)n,„ (6.20) 

and in particular this means Here we can again define some rate functions in the form 





= Nc 


mo 


^Nc + l 






rhrn 


= (iVc + l)K. + l) 




= A^c("m + l) 




(iVc + l)n,„ 



Rate7,7(K, x, ^ ^ d^Kz ^ (fk J cfk' (5(Atj2K(x) - lu)S(K - K2 - k + k') 



-Rnc 

X [l + /K(x)]/2(x)>V/(x,k)Wj(x,k') 



Rnc 



J d^K4 J d^k J (fk' (5(AwK4(x) - w)(5(K4 - K - k + k') 

X /k(x)[1 + /4(x)]>V/(x,k)Wj(x,k') (6.21) 



The term is then defined as 

a/(K,x) 



dt 



2 ij 



^ R'Atcij{K,x,ni{Nc) - nj{Nc))nimj (6.22) 



5. Three-quasiparticle terms 

,2 



Rate^+^(K,x,u;) = ^ rf^j^ ^ (fk' j d^k" J(AtJK(x) - tj)(5(K - k + k' + k') 

X [l + /K(x)]>V/(x,k)>Vj(x,k')>VK(x,k') (6.23) 

Rate^j|^(K,x,t^) = ^ d^j^ ^ d^k'^ d^k" (5(Aa;K(x) - tj)5(K - k + k' + k') 

X /K(x)W/(x,k)Wj(x,k')Wif(x,k') (6.24) 



The term is then defined as 

a/(K,x) 



dt 



{Ratc(+^ [K,x,1]a'(A^c) + ni{Nc) ~ nj{Nc)] nimjmK 
Rate[-]f [Vi.^MNc) + nANc) - nj{Nc)] m/njn/^}. (6.25) 



UK 
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The various rate in the growth, scattering, and three- terms 
are compatible with the same rates as seen by the condensate, 
i.e., equations (126,127,141,143,144) of QKIII, in the sense 
that the rate at which a process removes/adds a particle to the 
condensate-band matches exactly with the corresponding rate 
at which the process adds/removes a particle from the non- 
condensate-band. 



VII. CONCLUSION 

The principal results, in the sense that they can be applied, 
are the master equations of the time-independent local equi- 
librium approximation, Sect.|l|, and of the time-dependent 
local-equilibrium approximation of Sect.^ The first of these 
was derived in QKlll, but without the inclusion of mean-field 
effects. These are now included by means of the use of effec- 
tive potentials, which account for the majority of such effects, 
giving a description which is quite realistic. The use of this 
picture — with some technical simplifications — gives the ba- 
sis for the condensate growth methodology of [p^, and this 
seems to agree quite well with experiment. 

However, the assumption of a static non-condensate-band 
is not always reasonable, so the methodology of Sect.^ to 
get the time dependent local equilibrium approximation is 
necessary. Fortunately, the basic modification is essentially 
straightforward. The equations are modified so that the non- 
condensate band evolves according to a quantum Boltzmann 
equation, and the content of Sect.^ is how one includes from 
the point of view of the non-condensate-band equation of mo- 
tion the correct way of treating the transfer of atoms and en- 
ergy between it and th e condensate ba nd. The necessary mod- 
ifications are given in ( 5.12 , 6.22 , 6.25 ), are easy to understand, 
but have a complexity which arises from the mixture of parti- 
cle creation and destruction operators which defines the quasi- 
particle creation operators — this complexity is intrinsic to the 
problem and therefore unavoidable. 

The whole procedure has been developed in terms of the 
Bogoliubov picture of condensate-band atoms moving in the 
effective potential, which includes mean-field effects, and is 
essentially a kind of Hartree-Fock picture. Even though no 
explicit depiction of the scattering of quasiparticles by each 
other is given, this can be treated if required by one of the 
following options. 

1 . We can ignore these scattering effects entirely, when they 
are negligible compared to the scattering of quasiparticle by 
non-condensate-band atoms. This will certainly be the case if 
there is a very large proportion of thermal vapor 

2. Notice however that the condensate band Hamiltonian H(? 



as defined in (3.36 )is never itself explicitly approximated by 
the Bogoliubov form in the master equations. It is only the 
transition rates induced by collisions which have been ex- 
pressed using the Bogoliubov approximation. One could in 
principle solve the master equation using the full H^, and the 
nonlinear mixing induced by the nonlinearity must be equiva- 
lent to the collisions of quasiparticles, because of the presence 
of the randomization induced by the collisions with the vapor 



3. What seems simplest is a compromise — approximate 
by the Bogoliubov form taken to the next highest order, which 
gives rise to Belyayev terms. In the 1 /\/N expansion in which 
the Bogoliubov method is valid the next higher order terms 
must be quite small. This is possible because the the con- 
densate band consists only of a small fraction of the levels 
available, and the quartic terms which appear non-negligible 
in a conventional Bogoliubov [^] method are taken account 
of her by the use of the effective potentials. 

Thus the quantum kinetic methodology is complete. It pro- 
vides a logical combination of kinetics, as given by the quan- 
tum Boltzmann equation, with the field aspect, as given by 
the Bogoliubov picture. It is likely to be as accurate as any 
experiment will require, and has already been used to give a 
quantitative picture of condensate growth [ p"6HT^ ], consistent 
with experiment. 



APPENDIX A: CONSTRUCTION OF WAVELETS 



1. One dimensional systems 



In this case the wavelets can be expressed directly in terms 
of trap energy eigenfunctions. Consider the one dimensional 
trapping Hamiltonian 



Y^E^\n){n\ 



(Al) 
(A2) 



We relabel these states by dividing them into groups labeled 
by an index N = 0,1, 2,..., so that there are A^v states in 
each group, and we can now label the states for a fixed by 
a subsidiary index I. Thevalue of the parameter Aat will be 
fixed in due course, to give a description of maximum conve- 
nience. 

This means we write the states as \N,l)g, where I = 
1,2,3, An — 1, and now the energy of a given state is 
En{Nd), with 



JV-1 

n{N,l) = 1-1+ An'. 

N'=0 



(A3) 



Using this relabeling, we define the basic wavelet states as 



-7ri(Ajv + l) ^ 

\N, m) = J2 e'"""/^" \N, l)g, (A4) 



where m = 1, 2, 3, . . . , A 



1=1 



(A5) 
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N = 6 



N = 5 



N = 4 



N = 3 



N = 2 



N= 1 
i N = 



Fig. 2. Grouping of Trap levels. Wavelets are formed by linear com- 
binations of levels in each group. 

These states form an orthonormal set 



{N',m'\N,m) = Snn'^t, 



(A6) 



The actual wavelet functions are given by the position space 
representations of these states; 



WN,7n[x) = (.T|iV,m) 



(A7) 



{N,e\HT\N',6') 

= 5NM\EN6ee' + [{N,e\HT\N\e') - EnSso] } 

(All) 

It is now necessary to make some approximations in order to 
get some idea of the kinds of functions which will turn up. The 
major approximation is to assume the energy of each level is 
approximately a linear function of / for any N . This is exact 
in the case of a harmonic potential. 
Thus we write 



E[N,l) = E{N) + ei,{l) 
where now, we approximate 



(A12) 



(A13) 



and this equation defines ojn- After some manipulation we 
find that 

{N, 9\Ht\N', 6') = 5n.n' {EN5ee' - ihc^NhNiO' - 0)} . 

(A14) 



Here h^lx) is defined by 

d f sin(AAra:;/2) 

hN[x) - 



dx [AArsin(x/2) 



(A15) 



a. Phase space variable 9 

The best physical variable to use here is 
6 = 27rm/AAr. 
This variable is an angle with the range 

27r 2 X 27r 3 x 27r 



27r. 



Aat An Aat 
This range approaches a continuum as A^r ^ oo. 



(A8) 



(A9) 



c. Discrete delta-function approximations 

When X is an integer, the function 

sin (Ajvx/2) 



kN{x) 



Aat sin {x/2) 



(A16) 



is a Kronecker delta function S^fi. Correspondingly, when 
hiy{x) is used in a summation like 



(A17) 



b. Expression of the trapping Hamiltonian 

The average energy of the states | A^, I) in the energy band 
is 



-. Ajv 

Ej,^—Y,E{N,l) 



(AlO) 



We can now write the trapping Hamiltonian Ht as a part di- 
agonal in both N and 9, and a residue which is diagonal in N 
only; 



provided F{x') changes sufficiently smoothly as x' 
. . . ,x ~ 2,x ~ l,x,x + 1, X + 2, . . ., then we can write 



F{x-1)- F{x + 1) 



(A18) 



and thus h]s[{x) behaves like a "derivative" of a Kronecker 
delta function. If F{x) is smooth over a wide range of integer 
values of x, so that the ran ge of x can be regarded as a contin- 
uum, the finite difference ( A.18) becomes essentially the same 
as a derivative. This is the limit in which we shall use these 
functions. 
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d. Phase space behavior of the wavelets 

By using the WKB approximation, we can derive an ap- 
proximate form for the wavelet functions in the form 



1 



exp - 



^^^ sm{AN[e + e{x,EN)]/2} 
sm{[e + 0{x,EN)]/2} 

In this equation, the function 



(A19) 



(A20) 



is the classical momentum of the atom at point x, while the 
function 



9{x,En) = MujN 



dx' 



Pe„{x') 



(A21) 



is proportional to the classical time taken for the atom to reach 
the point x; indeed if the trap is exactly harmonic with trap 
frequency lu, then Tiujn ~ ff-^, and 9{x,En) is exactly the 
phase angle ujt of the oscillator at the time taken to reach the 
point X. 



The approximate form of the wavelet function ( A.19) shows 
that it is sharply peaked when 6 = —d{x,E]y), so that the 
wavelets represent a localization at the point where a classical 
particle would have phase equal to 9{x, En) 




Fig. 3. The phase space cells for a harmonic oscillator 



e. Harmonic trapping potential 

In the case of the harmonic oscillator potential V{x) = 
kx^ 12, we can label the wavelets by a momentum p{N,9) 
and position x{N, 9) given by 



r{N,9) = ^j2EN/k cos 9 
p{N,9) = \/2MEn sm9 



{All) 
(A23) 



and hence get a set of values of r and p which can be used 
to label the wavelets. This set is completely equivalent to the 
TV, 9 labelling. 



/ Size of the phase space cells 

The value of Aat is not yet fixed, and it should be cho- 
sen for the maximum convenience. The best choice is one 
which makes the phase space cells all have approximately the 
same dimensions. For large N, we will show that this can be 
achieved by the choice 



A 



AV-Bjv, 



(A24) 
(A25) 



With this choice, the uncertainties of the phase and energy 
variables are 



5En 
S9n 



hujAy/ En 
2tt 



so that, using ( i\22| , [\23 ), the uncertainties in r and p are 

27rV2 , ^, hujA 



(A26) 
(A27) 



Sr 



6p 



AVk 
27^^/2M 



|sin6l| 



2V2fc 



2V2 



(A28) 
(A29) 



Here Sr and Sp vary quite strongly with the phase 9, unless 
the coefficients of the sin 9 and cos 9 in each of these equa- 
tions are th e same, and this happens when A = 2\j2'!: jUuj, 
as in (|A25|). 



In this case the uncertainties become 

cos 6*1) 

bp = \/ 7r7i7;(| sin^l + | cos 6*1). 




(A30) 
(A31) 



with?; = {kMfl^. 

The coefficients of | sin 9\^\ cos 6*1 are essentially the vari- 
ances of r and p in the ground state of a harmonic oscillator — 
thus we find h/2 < 5r 5p < h, which means that the cells cor- 
respond basically to minimum uncertainty states. The "cells" 
are thus square when measured in these natural units, and the 
variation in 6r and 6p arises solely because their orientations 
depend on 9. 

The choice ( A24 ) of A^r naturally leads to the the fact that 
for large N 



En 



(A32) 
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2. One dimensional field Hamiltonian in terms of wavelets 



3. Three dimensional formulation 



The field Hamiltonian can be written in terms of the non- 
interacting and the collisional parts 



H — i^kin + H-Y + Hi. 

We can write the non-interacting part as 



(A33) 



NO 



Nee' 



(A34) 



Here Mn {6 — 0') is defined in terms of the matrix elements 
dAliljAl^as 



MN{e-e') = -ihujNhN{0-e'). (A35) 



a. Phase space description 

If we use the labelling in terms of p and r, the structure of 
the description is unaltered, but the appearance of the expres- 
sion looks different. Thus we can write 

-ffkin + -f^T = ^ ErpalpGrp + ^ M Irr' pp)alj,ar' p' ■ 

rp rr'pp' 

(A36) 



This expression conceals the nature of the wavelet basis, 
which is implied only by the fact that the phase s pace point s 
r,p represent wavelets through the definition (\22,A23), 
which arranges the points on discrete energy shells at posi- 
tions on these shells determined by the values of 6. The na- 
ture of the function M{rr'pp') means that it is nonzero only 
for points r,p and r' ,p' which are on the same energy shell. 
When acting on sufficiently smooth functions we can use the 
approximate relations of sect. |A 1 c| , and in particular (A18), 
F{r,p), and thus approximate 



M{rr'pp')F{r\p') w ih [ - kr-^ ] F{r,p). 



M dr 



dp 



(A37) 



Thus the Liouvillian flow arises for sufficiently smooth func- 
tions, and we can formally write 



— p d 3 



(A38) 



provided a,.p is interpreted as being evaluated in an average 
over a density operator for a state which satisfies the smooth- 
ness requirement. 



The case of a harmonic potential is particularly simple to 
treat, since the energies in the three different dimensions are 
separately conserved. We will use a notation in terms of these 
three energy components as follows. Each "cell" is specified 
by a set of six quantities (A^i, N2, N^, 61,62, 6^), from which 
we can form the two vector quantities r, Q, defined by 



k 



cos 6i, 



fiQi = J 2MEN^i sign (sin di) . 



(A39) 
(A40) 



The definition of Q ensures that the set r, Q do determine 
{Ni,N2, N^i, 61, 62,63) uniquely; the sign of Qi determines 
the sign of the i-component of the momentum, since for a 
given Ei and this can have two possible directions. In 
terms of these variables, we now use the notation wq (x, r) 
for the wavelet function, and the notation Aq ,. for the corre- 
sponding destruction operator, so that the field operator can be 
expressed as 



Q 



(A41) 
(A42) 



The commutation relations are of course 

= 5rr"5QQ' (A43) 

The commutation relations for the ijjQ (x) will be needed; they 



are 



[V'Q(x),^t^,(x')] = <5QQ.^^i;Jj(x,r)z«Q(x',r) (A44) 



3Q.Q' 



2 2 ; ^ ;gQ(x,x ). 



(A45) 



The approximation to get the second line comes from the 
WKB form, ^A^. 



a. General potentials 

The construction here is specific to a potentials which can 
be written as the sum of x, y and z dependent parts, of which 
the three-dimensional harmonic potential is the most useful 
example. Even though harmonic traps are most commonly 
used, the effect of the condensate on the other particles acts 
like an additional potential, and it will therefore be necessary 
for us to consider more general potentials. We will, without 
proof, assume that it is always possible to find an appropri- 
ate wavelet basis of the kind set up here. Although this may 
not be so, particularly in the case of non-integrable systems, 
it seems very likely that the results of our analysis will con- 
tinue to be true, we say this because our methodology really 
depends only on a local properties, and almost every potential 
we are likely to use is locally harmonic. 



24 



APPENDIX B: TIME DEPENDENT MEAN-FIELD TERMS 

We will show that, under the requirement that that the mo- 
mentum disrtibution function is smooth, we may use the ap- 
proximate expression 

(V'^(x)^.tj(x')) = e'K(Q-)-|.g(x,x')F(K(Q,u),u) 

-^(Vv,9(x,x')) • VkF(K(Q,u),u) 

(Bl) 

in which x = u + v and x' = u — v. We will do this explicitly 
for the case of no trapping potential as follows. In this case 
we consider 



(^j,(x)^i,(x')) 



(27r)3 



K+A 



^i(k-x-k'-x') 



(B2) 



We can now express the smoothness requirement as the con- 
dition that the mean (a^ak') can be expressed as a function 
-R(p,q) of p = (k + k')/2 and q = k — k'. The require- 
ment that the different bands in K be uncorrected, which is 
fundamental to the method, only makes physical sense if any 
correlation that does exist between k and k' is over a range 
in k that is very much less than A, the spacing between the 
bands. Thus i?(p, q) must be sharply peaked in q. On the 
other hand, the dependence of i?(p, q) on p represents the 
overall momentum distribution, which should be smooth over 
the range A between the bands. 

This means that we can change to the variables p, q in the 
integral, and to a good approximation let the limits of the q in- 
tegration go to infinity, while those of the p are K— A, K+A. 
We implement the smoothness in p requirement by using a 
Taylor expansion for i?(p, p), getting first order approxima- 
tion 



1 



(2^)3 



K+A 

3„ „-i(p-u+q-v) 



d'^q / d'^pe 

oo -/k-A 

x|i?(K,q)(p-K).VKi?(K,q) 



(B3) 



'^•^|.g(x,x')^^(K,u) -i(Vv.g(x,x')) • Vk^^(K,u) 



(B4) 



where in this case 

sin Ax{x — x') s'mAy{y — y') sinAziz — z') 



5(x,x') 



it{x — x') TT{y — y') 7r(z. — z') 



(B5) 



as in QKI. 

Using t he formula (B3), tjhe commutation relations 
(2.12,2.12) and the Gaussian factorization of the correlation 



functions, it is not difficult to derive the form of the correction 
(6.5). It should be noted that the smoothness assumption will 
not be valid at the edge of the non-condensate band, where 
a boundary layer will form. This results from the fact that 
the Liouvillian straming in the time-dependent potential is not 
parallel to the actual fixed boundary between the two bands. 
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